Renormalized Thermal Entropy 
in Field Theory 



Sergio Cacciatori a , Fabio Costa b ' c and Federico Piazza , 

a Dipartimento di Scienze Fisiche e Matematiche, 

Universita dell'Insubria, via Valleggio 11, 1-22100 Como 

b Institut fur Quantenoptik und Quanteninformation (IQOQI), Osterreichische Akademie der 
Wissenschaften, Boltzmanngasse 3, A-1090 Wien, Austria 
c Perimeter Institute for Theoretical Physics 
Waterloo, Ontario, N2L 2Y5, Canada 

Abstract 

Standard entropy calculations in quantum field theory, when applied to a subsystem of 
definite volume, exhibit area-dependent UV divergences that make a thermodynamic in- 
terpretation troublesome. In this paper we define a renormalized entropy which is related 
with the Newton- Wigner position operator. Accordingly, whenever we trace over a region 
of space, we trace away degrees of freedom that are localized according to Newton- Wigner 
localization but not in the usual sense. We consider a free scalar field in d+1 spacetime 
dimensions prepared in a thermal state and we show that our entropy is free of divergences 
and has a perfectly sound thermodynamic behavior. In the high temperature/big volume 
limit our results agree with the standard QFT calculations once the divergent contributions 
are subtracted from the latter. In the limit of low temperature/small volume the entropy 
goes to zero but with a different dependence on the temperature. 



1 Introduction 



Thermodynamics is a very powerful tool for describing complex physical systems. Beside its 
evident experimental success in the laboratory, thermodynamics is used in the everyday practice 
of cosmology, for instance when the conservation of entropy is applied to a comoving volume of the 
expanding Universe (see e.g. p.]); on a more speculative level, the quest for a theory of quantum 
gravity often makes use of thermodynamic and entropic arguments, since they are supposed to 
be somewhat independent of the details of the underlying dynamics. Especially in view of such 
adventurous applications, it is certainly worth understanding better and better the connection 
between thermodynamics and microphysics within those physical regimes which are best known 
and under control. The tendency to thermal equilibrium, that in classical statistical physics relies 
on ergodic or mixing hypotheses, in quantum mechanics seems to be naturally driven by the 
correlations that a subsystem inevitably develops with its environment. This general view, which 
has occasionally appeared in the literature at the level of common wisdom, is now being put on 
firmer ground e.g. in the interesting book [2] and in related ongoing works. During time evolution, 
the reduced quantum state of a subsystem tends in fact to approach that of a thermal Gibbs state 
Pthermai °c e~ l3H , H being the Hamiltonian operator of the subsystem. This happens, under some 
generic circumstances and on time scales which are thoroughly discussed in [2J, even if the initial 
state of the subsystem is very different from pthermai, e.g. in the case when the entire system is 
initially prepared in a product state (subsystem) ® | environment). 

Our most successful microscopic description of physical interactions, quantum field theory 
(QFT), faces some difficulties when asked to reproduce coarse-grained meaningful thermodynamic 
quantities. In particular, as first noted in [3], the UV- divergencies encountered in the calculation 
of entropy are of a relatively uncommon type. If a finite system is in a thermal state, its entropy 
can be calculated with standard methods giving a thermodynamically sound result (see e.g. [1]). 
For the reasons described above, however, it is also interesting to consider, instead of the entire 
system, a subsystem occupying a finite portion of the entire volume. In this case, the entropy 
exhibits a UV-divergent "vacuum" contribution proportional to the boundary of the subsystem. 

1.1 Thermal Entropy and UV Divergences 

To be more definite, consider a system S whose dynamics is described by a QFT Hamiltonian 
H and put it in a thermal state ptotai <X- e ■ Then consider a region of space P (P stands for 
"place") of definite volume inside S and call R the rest of the system. The state in P is obtained 
by tracing out the irrelevant degrees of freedom belonging to R, p = Tr^ptotai- Then calculate the 
Von Neumann entropy S = —Trp In p, which is the appropriate generalization of thermodynamical 
entropy for generic quantum states (see e.g. [SI [5]). Schematically, in four spacetime dimensions, 
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one finds 



S = S vac (A,A) + S therm (V,T). 



(1) 



Here A and V are the boundary area and the volume of P respectively, A is a UV cutoff and 
T = 1//3 is the temperature. The term S vac is the UV-divergent entanglement entropy of the 
vacuum (see, among others [H EJ [91 [10]), obtained with the same procedure in the limit of zero 
temperature, i.e. when ptotai = |0)(0|. The general form of S vac is [11] 



c being a regularization dependent number of order one, for more general boundaries see, e.g. 
[12] . The expansion (j2J) follows quite generally e.g. from the heat kernel methods used in [9]. The 
finite component Stherm^ T), on the other hand, is the meaningful thermodynamic quantity; for 
a massless field it typically scales as ~ VT 3 in the big volume/high temperature limit. A few 
comments are now in order. 

While the leading area-dependent vacuum divergence (j2j) can be checked in a variety of ways, 
bringing out the subleading finite term Stherm is not trivial and, to our knowledge, has been done 
explicitly only for conformal field theories. In 1+1 dimensions Calabrese and Cardy [13], by 
exploiting the analytic properties of the theory, found a structure of the type with a term S v3iC 
logarithmically divergent with the cutoff. By using insights from AdS / CFT correspondence, Ryu 
and Takayanagi have been able to extend the result to higher dimensional [2]. However, it would 
be surprising if thermal entropy did not have the structure (pQ) in general. The entanglement of 
the vacuum is in fact a UV effect and should be there also for generic finite energy states. The 
two separate terms in (JT]) are thus expected in any plausible QFT theory where the highest energy 
modes decouple from the low energy physics. Of course, if thermal entropy did not have the form 
([T]) and, say, some divergent term were also temperature dependent, then the intent of deriving 
meaningful coarse-grained quantities would be even more troublesome. 

One may object that the volume dependent term in ([1]) wins over the area-term in the ther- 
modynamic limit. In thermodynamics, however, volumes have to be big in comparison, say, to 
the typical distances between particles. On the other hand, by taking for instance a typical cos- 
mo logical setup (T 2 ~ Mpi x Hubble) and A ~ Mpi, it is easy to see that the thermodynamic 
term overcomes in ([T]) only for volumes much bigger than the Hubble scale! We expect thermody- 
namics to be applicable in much less extreme conditions. One may also object that, rather than 
the absolute value, only entropy differences are meaningful. Still, the area dependent term in ([1]) 

1 In theories with AdS gravity duals, the thermal entropy ((T|) is conjectured to be proportional to the area of a 
suitably defined hypersurface living in a AdS black hole bulk. The two contributions in eq. ([TJ) have then two clear 
and distinct origins: the portion of hypersurface closest to the boundary of AdS accounts for the divergent "vacuum" 
term; the finite part is ! due to the bending of the hypersurface deeper in the bulk, due to the presence of the 
central black hole. 
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spoils the attempts of a thermodynamic description for subsystems whose size is (adiabatically) 
changing in time. In this regard, once again, a comoving volume in an expanding Universe is 
perhaps the cleanest example. 

Due to the non trivial dependence of the divergence on area, the quantity ([I]) cannot be renor- 
malized by standard methods, i.e. by adding local counterterms to the Lagrangian. Moreover, 
since the result ([2]) has been carried out for free and conformal theories, we are bound to have 
divergences regardless of the asymptotic behavior of the couplings or the UV completion of the 
theory, as long as such a completion is still a field theory. Of course, as proposed in [3J, one can al- 
ways subtract the divergent terms. The latter, however, are not more "spurious" than the widely 
accepted entanglement of the vacuu In plain contradiction with the general view/common 
wisdom illustrated at the beginning of this introduction, equation (CQ) is just saying that the state 
of a generic subsystem in QFT is actually very far from being thermal! 

Since the Hamiltonian operator H is an integral of a local density, one may expect that a state 
of the form e~^ u would factorize over contiguous regions of space giving e~ ,3Hp ®e~ l3HR , where Hp 
and Hp are the integrals of Hamiltonian density extended only to subsystems P and R respectively. 
If this were the case, tracing over R would trivially give a thermal state in P. However, the energy 
H hides a relevant amount of non-extensiveness that doesn't allow this factorization. The "inside" 
and "outside" contribution, Hp and Hp, do not add up to the total Hamiltonian because of the 
UV-divergent contact term Hj coming from the gradients across the boundary between P and R. 
In QFT, because of the singular nature of the interaction Hi between P and its environment R, 
the general arguments of [2] are not applicable. 

1.2 Renormalizing the Entropy 

The above difficulties can be ascribed to an inconvenient choice of degrees of freedom. To see what 
that means, note that the system/region of space in question has two complementary descriptions 
[T6l \T7\ [18] . In compliance to common intuition, P is described in classical general relativity by a 
subset P - more specifically, a submanifold - of the points/events at a given time-like coordinate 
t. On the other hand, as a quantum subsystem, P is described by a Hilbert space Ttp, which 
is a factor in the tensor product decomposition (H = Hp <S> Hp) of the total Hilbert space of 

2 It is sometimes emphasized that the Von Neumann entropy of a subsystem has a different meaning depending 
on whether or not the entire system is in a pure state. Thus, one may be tempted to take the entanglement of 
the vacuum seriously [25] and just reject the divergent term S vac in {T]) as spurious. Consider, however, the entire 
Universe in a pure state; we can expect thermalization to occur over some region S after local equilibrium is reached 
i.e. ps ~ e~P H . Deep inside S one can eventually consider the subsystem P. The Von Neumann entropy of P 
is clearly both a thermodynamic entropy - because S is in a nearly thermal state - and an entanglement entropy 
- because the whole Universe is in a pure state. Note also that, starting from the entangled vacuum of QFT, 
one can always construct states which are less and less entangled; such states, although probably not generic nor 
thermodynamically very interesting, would end up having a negative entropy after the subtraction of S va . c . 



3 



the field theory under examination.! In this paper we show that if we associate the region of 
space P with a more suitable - although unconventional - set of quantum degrees of freedom Hp, 
the entropy is already " renormalized" and has a perfectly sound thermodynamic behaviour. For 
a massless free scalar in the large volume/high temperature limit we find S = c(d)VT d , where 
d + 1 is the dimension of spacetime and c(d) a numerical factor that we calculate explicitly (see 
eqs. (|67p ). In this limit our results are consistent with those obtained by [T3], [T3] once their 
divergent contribution is subtracted. However, there is no trace of any area-dependent term in 
our calculations and no infinities are encountered except, of the IR type, in the 1 + 1 dimensional 
massless case in the limit T — > 0. In the same limit, our results differ from the "finite piece" Stherm 
that is found in ([1]) by using the conventional approach, although entropies tend to zero in both 

cases. For T -> we find in fact S VT d \n(VT d ), while [Ml Q3] find S therm ~ {yT d ) i - d+1 ^ d . 

We comment on this in the conclusions. 

Before describing our calculation in section [2] it is perhaps worth spending few more words and 
see in which sense our anomalous renormalization procedure underlies a different "localization 
scheme" (see [18] for more details). The main point here is how to pick a bunch of (local) 
quantum degrees of freedom Hp out of a larger system. A tensor product structure (TPS) - or 
quantum partition - can be assigned by specifying the observables of the individual subsystems 
[T9l [20] . In a composite system Hp <8> Hp two sets of observables A j (P) and A k (R), separately 
defined in P and R respectively, commute by construction: 

[A j {P),A k {R)\ = for any j, k. (3) 

The point here is that such a trivial result can be applied the other way around [TH] : if, within the 
algebra of observables acting on H, we manage to isolate two commuting subalgebras A^(P) and 
A k (R), they induce a uniqucl bipartition H = Hp® Hp on the whole system. In the conventional 
calculation of entropy ([1]) it is implicitly assumed that the quantum degrees of freedom Hp of a 
region of space P at time t are those defined by the set of local relativistic fields <p{t, x e P) and 
their conjugate momenta 7r(t,x G P). In fact, thanks to the canonical commutation relations, the 
two subalgebras generated by and 7r with labels x inside and outside P satisfy ([3]), and therefore 
induce a TPS. Such a TPS is the conventional localization scheme in QFT. 

In order to renormalize the entropy, in this paper we use an alternative set of commuting 
operators - and their corresponding TPS - as a new rationale to isolate the quantum degrees of 
freedom of P. We consider a free scalar field in d + 1 dimensional Minkowski spacetime. The 

3 Note that partitioning a quantum system [TSJ, [TH] is more subtle than making the partition of a set. While 
a finite set admits a finite number of possible partitions, a finite number of quantum degrees of freedom can be 
divided in an infinite number of inequivalent ways. The simplest example is a four dimensional quantum system, 
H = C 4 , whose partitions into "two spins" , ~ C 2 <S> C 2 , are in one to one correspondence with the elements of the 
infinite group S77(4)/S*[/(2) 2 [TH]. 

4 Actually, only if the two subalgebras generate the entire algebra of operators on Ti [Jj|] 
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normal ordered Hamiltonian reads 



H 



I 



d kwkCL^ ak, 



(4) 



where Wk = V^ 2 + m 2 and operators satisfy the commutation relations [at, aw] = 0, [ak, = 
<5 3 (k— k'). Instead of the relativistic fields and their conjugate momenta we introduce the "Newton- 
Wigner" fields a(x) just as the Fourier transforms of a^. 



The above defined operator a^(x) is directly related to the Newton- Wigner (NW) position operator 
[2T] in that, acting on the vacuum, it produces an eigenvector of eigenvalue x. Note that the 
relativistic invariant measure l/y/2wk is absent from the integrand and therefore those operators 
are not relativistically invariant. This amounts to the fact that a particle perfectly NW- localized 
according to some observer is instead "spread" when described by a boosted one [22] . On the other 
hand, the dynamics is still relativistically invariant because we are not changing the Hamiltonian 
of the free scalar nor the other generators of the Poincare group. What is usually considered a 
drawback of the NW approach doesn't worry us too much, particularly here since thermal states 
with T > break Lorentz invariance. The obvious reference frame in this setup is the one with 
four- velocity parallel to the expectation value of the momentum operator, . We refer to the 
extensive literature for more technical details (e.g. [21], [22], [23]) and philosophical implications [22] 
of NW operators. A very introductory comparison between the two localization schemes is found 
in [E]. 

2 Sketches of the Calculation and Main Results 

Since the entire problem is stationary, we can forget about the time coordinate. At some given 
time we cut a. d + 1 dimensional Minkowski space into two connected regions: P, of finite volume, 
and R, i.e. P U R = M 3 . Our results are independent of the shape of P. We distinguish spatial 
coordinates belonging to different regions using labels p, p', Pj ■ ■ ■ for points inside P, r, r', r j . . . 
for those in R and x, y . . . for generic points in ~R d . One of the basic properties of NW localization 
is that the vacuum of the theory is a product state, i.e. |0) = \0p) ® \0r) [18]. Moreover, starting 
from the vacuum, we can repeatedly apply a^(x e P) and a^(x G R) and generate two independent 
Fock spaces: 



This Fock space decomposition of the regions is the distinctive feature of NW localization and 
allows an intuitive representation of particles localized in space. This is impossible in the standard 





Hp = C © P x © P 2 © . . . © P n © . . . 



Hr = C © i?i © R 2 © . . . © R„, © . . . . 



(6) 
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localization, as the Reeh-Schlieder theorem forbids the existence of particle states of finite energy. 
The most striking consequence is that the vacuum is here a product state, while it is entangled 
in the standard scheme. In each Fock sub-space of given particle number we choose the obvious 
basis 

P n -> |pi . . .p n ) = — L=a f (pi) . . .a f (p n )|0) R n -> |ri...r ft ) = -=a + (ri)...a t (r,i)|0). (7) 



When we calculate traces we therefore sum on a basis in a Fock space of given particle number 
and then sum over all Fock subspaces. For example, if we restrict to block diagonal states (as are 
all those appearing here), we have 

Ti R ■ = (0 R \ ■ \0 R ) + [ d d r{r\ ■ |r) + / d d rd d r'{vr'\ ■ |rr') + . . . (8) 
Jr Jr 

(here and in the following we write P, R, for P, R, understanding the identification of regions 
with subsystems via the NW scheme). 

We put the entire system in a (non normalized) thermal state, p to tai = e~^ H , where H is defined 
in We call p the reduced state in P: p = Trpp to tai and we calculate Von Neumann entropy 
by means of the formula [9] 

S = -Tr P (plnp) = {--^- + 1^ InTrpp™ (9) 

which allows us to lose track of normalization factors. 

It's easy to check that p (as well as ptotai) is block diagonal on the particle number subspaces 
of its Fock space. Its matrix elements on each subspace are expressible in terms of the crucial two 
point function defined in the one-particle sector 

A-(p,p<) S M^), (10) 

where Q = (0p| p |0p) is the vacuum-vacuum matrix element. The function K is itself an infinite 
series obtained by tracing over R (see eq. (|24p below) and can be nicely written (see eq. f[iZ6T) ) 
in the diagrammatic formalism to be introduced in the next section. Since we are dealing with a 
free theory, the generic matrix element (P1P2 • • • Pn|p|p'iP2 • • • Pn) i s expressed as a combination 
of products of Ks in equation (|29|) . 

Since p is an operator acting on P, when calculating matrix elements of p n , integrations over 
the variables p have to be carried out in each Fock subspace. A final integration over the same 
variables has to be done in order to obtain Tr pp n . The terms in the corresponding series rearrange 
(see eq. fl33l) ) and this number can be written in closed form, again, in terms of the two point 
function K as 

Tr P p n = n n exp [ V -TrK j - n ] , (11) 
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where, by definition, 

TrK m = J^d d p 1 d d p 2 ...d d p m K(p 1 ,p 2 )K(p 2 ,p 3 )...K(p m ,p l ). (12) 

Note that the normalization factor Q n in (ITT]) drops when used in formula (Q . Note also that each 
factor K in the above integral is itself a series of integrals over the r variables. A consistent part 
of this work is finally devoted to evaluate the above quantity in the high temperature, V/ [3 d — ► oo, 
and small temperature, V/(3 d — > 0, limits. Two distinct behaviors of Tr K n as a function of n 
follow. 

In the high temperature limit the leading term of the series giving Tr K n is the one containing 
only integrations over the p variables: the other terms converge and give a subleading contribution. 
The actual proof of that is rather involved and part of it appears in the Appendix. The only 
integrations left are those inside P. Those are all of the same order in VT d , although they scale 
as l/n d , where n is the number of integration variables. The corresponding behavior of Tri^ n , 
when used in (fTTT) and (J9]), gives an extensive entropy. In this limit, apart from numerical factors 
to be found in eqs. (|66j) and (1671) . one finds in fact {l/VT d — > 0) 

VT d 

TrK n ~ —r + 0(l/VT d ) , S ~ VT d + 0(l/VT d ) . (13) 

n d 

In the low temperature limit the external integrals are no longer negligible and they have to 
be summed up. The corresponding series can be explicitly calculated at leading order in VT d . On 
the other hand, terms with a higher number of internal integrations are subleading and this gives 
(VT d -> 0) 

TrK n ~ (VT d ) n + 0(VT d ) n+1 , S ~ -VT d lnVT d + 0{VT d ). (14) 

The numerical factors are found in eqs. (1751) and (176|) below. In section 14.3.21 we will also consider 
massive fields; in the low temperature limit they behave like (1T4"|) except that the quantity VT d is 
each time suppressed by the factor e~ m l T (see eq. ( 183"]) below). 

3 Formalism and Diagrammatic 

After the general setting described above (basically, eqs. (JIJ-(J7J)) we give here more details of 
the calculation. First of all, we need to calculate the matrix elements of ptotai — e~^ H in position 
space; using the basis vectors (J7J and going to Fourier space, we have 

(xi . . . x n | e~ pH | x ; . . . x^> = l _ - x -(i)) ■ • ■ ^ n - x -w)- ( 15 ) 
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Here, of the total n points, m is the number of points inside P, S n denotes the group of permu- 
tations over n elements and 

//j(x-x') = (xle^lx 7 ) = -^p J Ae lk - (x - x,) -^ . (16) 
Explicit expressions of the two point function / and its massless limit follow for d — 1 

W = ^fo *^v^) - ^f^pj + £V> (IT) 

and d = 3: 

« X > = ^(^) - ^Jj^pp + 0(m>). (18) 

Here i£j are the modified Bessel functions of the second kind. Note that ffl6|) is not the usual 
QFT thermal correlator and, as such, it is not periodic in (3. This reflects the fact that we are 
not working in the usual thermic representation, where traces are taken by functional integration 
over a compactified Euclidean manifold. 

A crucial property of the two point function Ip that follows straightforwardly from its expres- 
sion (Tl6l) in Fourier space is 

j ^/ /3 (x-z)/ 7 (z-y) = l3 +7 (x-y). (19) 

By iteration we also have 

i?(x-y) = / n/J (x-y), (20) 

where the n th power of Jg has been implicitly defined in an obvious way. 

Matrix elements with a different number of particles on the two sides vanish, because in 
transformation (J5J) the number of creation and annihilation operators is preserved; our matrix 
Ptotai is thus block diagonal in the subspaces of given particle number. The same property is 
retained by the reduced density matrix p with respect to the local Fock space Tip, so that we only 
need to calculate the matrix elements (pi . . . p n |p|pi • • • p' n )- 

First, we define Q as the matrix element of p on the local vacuum in P: 

n = (Op | | Op) = 

= ((0p|®(0p|)p(|0p)®|0p» + / tfV«r|® (0 P |)e-^(|0 P ) ® |r))+ 

[ d d r 1 d d r 2 ((v 1 r 2 \®(0 P \)e- /3H (\0) P ®\v 1 r 2 ) P )+ ... (21) 

J RxR 



s 



In terms of the two point function Ja(x — x') we have 



n 



oo _ „ n 

X)zr / d d n...d d r n U 1 ^-^ = 

HI I Dn 

n=0 JH o-e5„ j=l 

\ I d d ri d d r 2 [I p ( ri - r 1 )I^(r 2 - r 2 ) + Ip{r x - r 2 )J^(n - r 2 )] 

^ JRxR 



d d r Ip(r-r) 

* JRxR 

+ ... (22) 

We can write this kind of expressions in a diagrammatic form; in this way, the vacuum expectation 
value is given by the sum of all the "bubble diagrams" : 

\ 

f \ I / \ / \ / \ 

Q = 1 









+ . . . 



(23) 



/ 



Here and in the following empty circles o are points in R and full circles • points in P, lines are 
the two point function / and two lines getting at the same circle imply integration. The n th term 
of the series (1221) is obtained diagrammatically by taking n empty circles and connecting them 
with each other in all possible ways such that each circle is reached by two lines. 
The matrix element of p living in the one particle sector is a two point function: 

(pipip') = - p') + / d d r - p')i>(r - r) + J^(p - r)/^(r - p')] + . . . (24) 

JR 

diagrammatically, 

/ 



(pIpIp') 
/ 

l 

2 - 

V 



V 




\ 



-e- 



/ 



00 



o 



-e- 



+ 




+ 



(25) 



Each term, weighted by a factor 1/n!, consists of all possible ways that the two external lines 
with the full circles can connect each other through n empty circles. Note that the "vacuum 
contribution" factorizes out, leaving 



/ 



\ 



(pIpIp') = n 



+ 



-e- 



+ 



-e e- 



+ . . . 



(26) 
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We call the two point function inside the parenthesis K: 



(p\p\p') = ClK(p,p'). (27) 

When we consider the matrix elements in the n-particles sector, bubble diagrams again factorize 
out, leaving us with an expression of the form 

(P1P2 • • • Pn|p|p'iP2 • • • Pn) = ~ r ~ points connected diagrams) , (28) 

where now we are summing over all the possible diagrams that connect n points on the left to n 
points on the right. As we can have only two lines starting from each internal point, a diagram is 
composed by "paths" , each of which connects one point on the left to one on the right (two points 
on the same side cannot be connected). This means that, if we select a pair of points (pj, p^), we 
can factorize an expression equal to the sum of connected diagrams in fl26|) . which gives the two 
point function K(pi, p'-). To obtain all the diagrams, we have to consider all the possible pairs; 
the result is that we can write all the matrix elements in terms of functions K(pi, p'.): 

Cl n 

( Pl p 2 . . . Pn|p|p'lP 2 ■■■Pn) = — { Yl II K (Pi> P -0)) ' ( 29 ) 

' o-GSn 3=1 

Note that this expression has the same structure as f|T5|) . with /g(x — x') replaced by K(p,p'); 
the reason is that, using NW, local regions have the same Fock structure as the global space. 

We define the powers of the two point function K by multiplying K as a one-particle operator 
acting inside P. Accordingly, we define the trace of some power of K as 

TrK m = f dpdp' ...dp (m - l) K(p,p')K(p',p")...K(p^ m ~ 1 \p). (30) 

J pm 

These traces are what we need to calculate Trp p n which, in turn, allows us to find Von Neumann 
entropy by means of eq. [9] that can be applied to p regardless of its normalization. 

Consider first the case n = 2; as p is block diagonal in the fixed number of particles subspaces, 
so is p 2 , and we can write, for the generic matrix element in the m-particles subspace 

(PlP2..-Pm|p 2 |pip 2 ---Pm) 

= / rf d gi...rf d g m (pip 2 ...p m |p|qi...q m )(qi...q m |p|pip 2 ...p^) 

J pm 



ci- 2 



m r 

Q - 2 
ml 



d% . . . d d q m II K (P*' C l-W)^( c b-' Pa'(j)) 

<7,a'eSmij=l 

/ d d qi ...d d q m V \\K{p j ^ a(j) )K{(i (Tl{j) ,p' j ) 

j pm 
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o 2 m r 



a 2 



m! 



(31) 



CT eS m .7 = 1 

Iterating this procedure we obtain the expression for p n : 

(PlP2-.-PmP PlP 2 ---Pm) = — f 

m! 

Finally, the trace is given by the sum of all the contributions of the m-particles matrix elements: 



£11** 



(32) 



+oo 1 „ m 

m=0 m ' <rG5 m 7=1 



= fi n det(l -7T 1 )- 1 
(see for example [2E], page 187, formula (4-86)), and then 



(33) 



Tr P p n = Q n exp 
Inserting this expression in we find 

S = 



oo \ 

E -Tr if jn 



(34) 



d 

dn 



\ 1 
1 J E- Tr ^ 

' 7=1 ^ 



(35) 



n=l 



The quantities we need to calculate are Tr K n ; we can give for these a diagrammatic expansion, 
like in fl26l) . For n = 1, we only have to "close" each diagram, so to match together the full circles 
at the ends of each factor, and we get 






Tr/V V J ■ { J + V ;+..-. (36) 
Operator K 2 is obtained multiplying term by term two copies of the expression for K: 



( 



\ 



( 



m • ' • o • 



+ . . . 



m • ' • o • 



/ 
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— • • • + 2 « • e » + « e • e • + • • • ; (37) 

when we take the trace, again we have to match the extremes of each diagram, so that we end 
up with closed loops with two full circles each and an arbitrary number of empty circles. The 
analogue expression for TrK n is a straightforward generalization: it contains loops with n full 
circles and arbitrary empty circles. The explicit formula is 



d d Pl ...d d Pn / rfV i ■ ■ ■ rfr'frdtr'l cfVi n) • • • d d rf] 

Ip( Pl - vDW, -rl)... I p {vl - p 2 )J^(p 2 - r?) . . . I p {vl - Pl ) . (38) 



4 Explicit evaluations of entropy 

We consider separately the two situations (3 d <C V (high temperature) and (3 d ^> V (low temper- 
ature), where V is the volume of the region P under consideration. We will mainly consider a 
massless field, but we will also consider a finite mass m in the low temperature limit. 

The high temperature limit, V/f3 d — > oo, is also the limit of very large volume, so we may 
expect to find, at leading order, the entropy that we would find by considering as our subsystem 
the whole space. Although this turns out to be the case, the actual proof is pretty involved and 
is carried out in the subsection 14.21 



4.1 The Whole Space 

As a first check of our formalism we calculate the entropy of a system (not a subsystem) in a 
thermal state. This can be done by standard methods, i.e., by calculating the partition function 
Z = Tie- pH that, in our formalism, using eq. (fl5j) . reads 



+oo „ 1 n 1 

z = / d d Xl . . . d d x n - n w** - x -(i)) = ex p(E - TrJ ^) • 

n=0 J ' u£S n j=l j=l J 



(39) 



We basically used the same derivation as for (1341) . In fact, in this case, K = Ip since there is 
nothing to integrate over outside the system. The operator powers Ip are obtained by integration 
over all space. For this purpose, we use eq. (I20j) and find 



Tr/S = / d d xI n p(-K - x) = VI nP (0) . (40) 
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On the other hand, for m = 0, we have 

tjq) - [ ddk r -(>M - Qd r°° dkk d - i e-' 3k - Qd ^~^ ! (ad 

where = 27r d / 2 /T(d/2) is the d-dimensional solid angle. From (139|) we then find 

where £ is the Riemann zeta function. If we calculate the entropy by using 

S = (-f3^- + l)\nZ (43) 



d(3 

we find the leading order (1671) of the more general result in the high temperature limit. 
4.2 High temperature limit 

We consider the limit (3 d /V — > in the massless case. Let us first consider the one dimensional 
case where P is the the interval (-L, L). 

4.2.1 One dimensional case 

We need to compute TrK n . For n = 1 it is given by (136|) . By eq. (IT7|) the first term is simply 

J-L ^(P 2 + (X ~ X) 2 ) 7T/3 

where the term in parenthesis after the diagram, here and in what follows, specifies the suffix (3 
of the corresponding Ip functions running in the loop. For the second term, we have 







(P,0)= V J (2/3) " V J & ft 



1 dx [ dy——4- — ; f — "TtT • (45) 



7T2/3 i_ L J_ L »n(^ + (x- yf) vr(/5 2 + (y - x) 2 ) ' 
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In the first term of the latter equality the property (|20|) has been used. Now 

L f L (3 (3 

dx / ay- 



is J-L 



tt((3 2 + (x- y) 2 ) tt(/5 2 + (y- x) 2 ) 
(3 d f 2L J f 2L J 1 

dx I d v , — — 



2ir 2 d(3j J *(3 2 + {x-y)< 

(3 d f 2L , 1 x 1 f 2L , d f x 

/ dx— arctan — = — — ; - / ax— x arctan — 

vH d/3 Jo (3 (3 (3tt 2 J dx \ (3 

2L 2L L 

arctan — = — + 0(1) , (46) 



(3n 2 (3 n(3 




thus 



0(1) . (47) 

Thus it seems that terms containing external integrations give finite contributions. Before checking 
this to the next order, let us note that 

L f L m m 1 ---m n+1 /r+ 2 

dxn ■ ■ ■ I dx r . 



7T 



L J-L 



2L 



[f3 2 ml + (x - sx) 2 ] ■ • • [I3 2 m 2 n+1 + (x n+1 - x ) 2 ] 



n+2 



POO 


poo 




/ dx\ ■ ■ 


■ / dx n+ i 




Jo 


Jo 






momi • 


■ ■ m n+ i 



[m 2 + x\] \m\ + (xj - x 2 ) 2 } ■ ■ ■ [m 2 + (x n - x n+1 ) 2 } [m 2 n+1 + x 2 n+l ] J 
+ cyclic terms + . . . (48) 

Here cyclic means with respect to the dependence on rrii and the ellipses mean terms of the next 
order in /3/L (when (3/L goes to zero). This can be easily verified as follows: 
Let us call it U(mo, . . . , m n+ \). One must show that the limit 

lim yU(m , ■ ■ ■ ,m n+1 ) (49) 

p — >0 l_j 

exists and is given by the above expression. To this end, it is convenient first to rescale all 
coordinates by (3, so that 

j-U(m , . . . ,m n+ i) = 

- ft dr ■ ■ ■ ft dr m omi ---m n+1 

~ * n+2 L 7-| X ° 7_| ^ [m§ + (x - Xl ) 2 } ■ ■ ■ [m 2 n+1 + (x n+1 - x ) 2 } ' [M) 
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Setting z = L//3, it suffices to compute the limit z — > oo by means of the de l'Hospital rule, to 
get the desired result. Indeed the de l'Hospital rule says that we must look at the limit 

1 



lim 



k .k 

f 

""" Jl n+1 [ml + (xi - x 2 ) 2 } ■ ■ ■ [ml + (x n - x n+1 ) 2 } 



m^rni ■ ■ ■ m n+1 



1 



+ 



1 



1 



m l + (§ " x i) m l+i + (§ - x n+i) ™o + (§ + »i) rn 2 n+l + (f + x n+1 X 
+ cyclic, 

where we used the obvious relation 

d_ 

dz 

The two factors in the square brackets give the same contribution^!, so that after the shift 
Xi + L//3, for all the Xi, our limit becomes 

2 



f(x)dx = f{z) + f(-z). 



(51) 
(52) 



lim 



z— >oo 77 



n+2 



2L 

8 



dx\ ■ ■ 



dx 



m m 1 ■ ■ ■ m n+ i 



n+1 r™2 



ml + (x 1 - x 2 ) 2 ] ■••K + (x n - x n+1 ) 2 ] 



1 



m o + x\ m? nJrl + 



cyclic. 



(53) 



It remains to show that indeed the integrals on the r.h.s. of (|48l) converge. This can be done by 
introducing the new variables tj such that 

ti = Xi — Xi + i , i — 1, . . . , n , t n+ i = x n+ \ . (54) 

Then, after rewriting the integral, one sees that the integrand is dominated by Y\a=q ™i t 2 • 

i ' i 

Now, let us continue our analysis and consider the term 



(/3, (3,(3)= K^J (3/3) - Q (2/3, (/3, 2(3) + Q (/3, (3, (3). 
Using our general formula, find 





(55) 




(3/3) 



2L 

3^/3 



+ 



Q(2/?,/?) =0 



(A 2/3) 



8L 



1 



/3tt 2 7 (l + x 2 )(4 + x 2 ) 



+ . . . 



_8L_ 

3^/3 



1 + x 2 A + x 2 



+ ... 



2L 

3^/3 



+ 



3 the second one is obtained by the first one by changing sign to all the integration variables 
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0„ T /*00 poo 
(/?,/?,/?) = — / dx I dy- 



0^ J 



1 + x 2 )(l + y 2 )[l + (x - yf 



+ ... 



2L 

3n0 



+ ... , 



so that 



lim 



0^) (0,0,0) 



= . 



(56) 



(57) 



This is true for any power of K. For example 



(58) 



and using the above results we see that 



Q (0, 0, 0) = Q (0, 0,0)= Q (0, 2 0) ~ O (A P> ®> 



(59) 



and then the L/0 terms drop out. 

Indeed, one can prove that any loop integral containing at least a white ball insertion converges 
in the high energy limit. A complete proof is very tedious, but tracks of a proof can be found in 
the Appendix. 
Thus we get 



TriT 



Q 



(ra; 0, . . . , 0) + . . . 



(60) 



that is the leading contribution is given by the loop with exactly n black ball insertions. Now, 
being cyclic terms all coincident, 



Q(n; /?,...,/?) 
2nL r °° 



0* n Jo dXl "'Jo dXn - 1 (l + xl)(l + x i n _ 1 )[l + (x 1 -x^] 



+ 



2nL 

07l n 



U{n) + 



As we will see soon, the integral can be computed to give 



U(n) 



IT 



n-1 



[1 + (x n - 2 - Xn-i) 2 } 

(61) 
(62) 
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so that 



Indeed, we can show this as follows. Using the results shown in the Appendix, we can add a loop 
with a white insertion and n black insertions, without changing the divergent term, so that 

Q(n;/3,...,/3)~ Q (n; ft ...,/?)+ Q (n; A /?) 

= Q(n-l;2/3,£...,/3) , 
where ~ means equal up to convergent terms. Next we can add the term 

Qi(n-l;2/3,/3,...,/3) 

to obtain 

Q(n;/3,...,/3)~ Qj (n - 2; 3/3, ft . . . ,(3) . 
Proceeding in this way, we get 




(n;/?,...,/?)- \^J{n(3) , (64) 
which gives (|63|) . 

Inserting (163]) in (135]) . we get for the entropy 

s =w + -- (65) 

Note that the above result is finite and extensive; no subtraction needs to be made. 
4.2.2 The d dimensional case 

The computations detailed out for the one dimensional case can be extended to any dimension 
d > 1. Indeed, using ( fl6l) . one finds again that the loops containing only black insertions dominate 
so that (!60j) and (|64l) continue to be true. On the other hand 

(n/3) = . V Q .. + . . . . (66) 
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Thus we can easily compute the entropy, whose result is 

* = kz&m ld+mi+1)+ --- (67) 

It can be useful to separate the odd dimensional cases from the even dimensional cases. Using the 
duplication formula for the Gamma function one finds 

V (2k — 1 V 



V_ (k_+ 1)! 

f3 2k + 1 ~K 

For example, for 3 spatial dimensions we find 



■W. = o W \ M 2 <(2k + 2) + . . . . (K 



4.3 Low temperature limits 

In this section we consider two different low temperature limits. In the first case we take v — > 0, 
with 

V 

v = 



(2vr/3) d ' 

for fixed values of the product (3m between the inverse temperature and the mass of the field. 
This can be thought as the small volume limit at fixed values of mass and temperature. However, 
we can also interpret it as the low temperature limit of a massless field; thus the mass plays the 
role of a IR regulator so that successively the quantity (3m must be set to zero. 

In the second case, we consider the limit (3 — > oo for fixed values of the volume V and the 
mass m, that is the low temperature limit of a massive field. 



4.3.1 Small volumes and massless limit 

Looking at v — > as a small volume limit, the result may depend on the shape of P and the way it 
"shrinks" . To avoid this problem, we will suppose for P to shrink down isotropically. This means 
that for any m dimensional section of P, having volume V m , the quantity V m /f3 m must tend to 
zero when v — > 0. Note that such condition is automatically satisfied in the low temperature limit 
interpretation, if P is contained in a compact region. For integrals over regions shrinking down 
isotropically, we can then use the approximationj 

J p cfW(p) = (2vr)V(0) + CV). (70) 

6 this is the form of the spatial integrals after a substitution p — ■> p//3; the /3 d extra- factors are reabsorbed by 
rescaling the momentum variables as k — > /3k. 
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To estimate TrK n we must approximate all loop integrals. The simplest one, is the tadpole 
integral: 



Q = J dV^(p-p) = VI p (0) =v J d d ke~^ 2+ ^ 2 



(71) 



Using (I70p . we see that each integration over P can be simplified giving a factor v. For example, 
for the two-points loop we get 

= J d d pd d qla(p - q)Ip(q - p) 

= v 2 (J dtke-^+V™) 2 ^ + 0(v 3 ) . (72) 

Note that the last integral depends on the temperature only via the fixed product (3m, so that the 
resulting expression is of order v 2 . In the same way we conclude that the loop with n full circles 
is of order v n . Diagrams with empty circles can be computed using the identity 



d d r = I d d x- I d d p . (73) 

R 



The integrals over the whole space can be reduced using ffl9l) - fr20l) , whereas each integral over P 
gives a contribution proportional to u, which is therefore a subleading term with respect to the 
first one. For example, for the "mixed" two-points loop we get 

= / A / d H(p-*-p) - = 

J d d pl 2l3 (0) + 0(v 2 ) = £°° dkk d - l e^ k2+ ^ 2 + 0{y 2 ) . 

In conclusion, to estimate the generic loop integral can simply remove each empty circle by means 
of the substitution f R d d rIp(p — r)/s(r — q) — ► Ig + s(p — q), so that the leading order of each 
diagram is determined by the number of its full balls. 
By applying this procedure to all terms in fl36|) . we find 

TrK = V J ij9 (0) + 0[v 2 ) = y fi d V - / dkk^e-V^+UP™) 2 + (v 2 ) = uQ d C d + 0[v 2 ) , 
i=i 3=1 3 Jo 

(74) 
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where C d is a multiplicative factor that doesn't depend on v. 

The diagrammatic expansion of r YiK n contains n full circles at any loop, so that for the leading 
term we get 

+00 „ 

TrK™ = dd P' ■ ■ ■ dd PnI 3 AVi - P2) ■ ■ ■ hAVn - Pi) + 0{v n+l ) = 

ji...j n =0 J 
n / +00 \ 

II E ^(°) + °^ n+1 ) = (ySl d C d ) n + . (75) 

1=1 \j l= J 

Plugging this result into ffTTj) and then into (jU]), we find for the entropy 

S = -vtt d C d \nu + 0{y) . (76) 

Note that, for d > 1, is a finite number. Indeed, it is a function of /?m defined by the series 

+°° 1 p+00 

c d = Y,- d / dkk^e-y/^+ue™) 2 , (77) 

whose terms are positive and bounded by the terms of the converging series C d (0). Thus 

C d (J3m) < 0,(0) = J2~1 dkk d ~ l e- k = {d-l)\({d) . 

j=1 J Jo 

An IR divergence appears in the onedimensional massless case. Obviously it could be cured by an 
IR cutoff that limits the integrations to the external region R. However, a natural regularization 
is provided by the mass term. Indeed, we can fix (3m at arbitrarily small but positive values, so 
that the general term in the series (I77p satisfies 

+00 rjPm r+00 



dke -V^TU^ = / dke - y/V + Wmy + / dke -^ k 2 + {jPm? 

JO JjPm 



-00 



< e-^ m j/3m + / dke~ k = e"^ m (j/3m + 1) , (78) 

Jj/3m 



ensuring convergence. For d > 1, no divergences occur and in the massless case the entropy (176 
reads 



S XlnX , 

X = L d 2^T(^)C(d) (d>l,/?ro = 0). (79) 



Again, the above result is finite and no subtractions need to be made. 
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4.3.2 Low temperatures for a massive field 

Let us first consider the one dimensional case. For f3 — > oo, the two point function Ip(x) behave 
as 

/,(*) ~ {^f m ^ ■ (80) 
When x falls inside P, the whole expression can be approximated by a constant: 



hip) ~ W) - J^' mP > 

and each integration over P contributes with a term proportional to the volume; thus we get 



m e -mP 



Q = J dplp(p -p) = 2L I p (Q) ~ 2L 
Q(n points)- J ^dpi... J ^d Vn (Ip{$)) n ~ 



It follows that a loop with n + 1 integrations over P is subleading with respect to one with n 
integrations and we can use again the usual trick to substitute integrations over R with integrations 
over the whole space, finding 



n=l 



Now, each addend l n p(0) is suppressed by a factor (e m / 3 ) n ; so that only the first term in the sum 
is relevant when (3 — » oo. The same argument can be applied to Trf^ n , giving 



TiK n ~ ( 2L.I 7 ^—e~ m(B ) , (82) 



2np 

from which we get the following expression for the entropy 

S XlnX , 

x = 2L 0p e ^- (83) 

These analysis can be easily extended to any dimension. Indeed, the two point function Ipijp) is 
still suppressed by a factor e _m/3 for p e P, so that all computation work essentially in the same 
way, providing for the entropy an expression of the form fl83|) . where the exponential multiplies a 
function of (3 with an at most polynomial growth. We will not enter into more details. 
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5 Conclusions 



In this paper we have considered the problem of computing the entropy of a subsystem confined 
in a finite volume region of a quantum field theory system. As hinted in the introduction, the 
subsystem under consideration is taken to be deep inside a much larger system which has reached 
thermal equilibrium and is therefore described by a thermal Gibbs state. The entropy of the 
subsystem results to be divergent in a non standard but curable way. The divergent part is 
the entanglement entropy of the vacuum and is a function of a cut-off and of the area of the 
surface bounding the region. Here we have highlighted the general issue of assigning appropriate 
quantum degrees of freedom to the considered region. Indeed, the inside/outside separation is 
traditionally realized through the usual localization prescription, which attributes to that region 
the local relativistic fields therein defined, together with their conjugated momenta. This in fact 
leads to a tensor product separation (TPS) of the whole Hilbert space, but it must face the above 
described problem of the infinities in the calculation of the entropy. 

In this paper we have adopted the Newton- Wigner (NW) localization scheme to define the 
TPS of the quantum fields associated to classical spatial separation. We made use of the creation 
and annihilation operators which act on the vacuum generating the NW position eigenvectors. At 
first, one could feel disturbed by the fact that these are non relativistically invariant operators. 
However, this does not affect the relativistic characterization of the dynamics. It simply means that 
localization is not observer independent. After all, a place/region of space is always to be defined 
on some chosen t = constant hypersurface. Moreover, Lorentz invariance is naturally broken by 
thermalization. In the NW prescription the vacuum is a product state, so that the corresponding 
divergent, area-dependent contribution to the entropy never appears in the calculations. We 
argue that the NW prescription is the right one in order to treat the problem of coarse-graining 
a thermalized microscopic system. The first point we bring to support our thesis is the finiteness 
of the results we obtained using such a prescription. We worked with a Klein-Gordon scalar 
field on a flat Minkowski spacetime in arbitrary dimensions, and computed the Von Neumann 
entropy for a subsystem confined in a region of volume V at a finite temperature T. The NW 
prescription automatically regularizes the ultraviolet divergences and directly gives finite results. 
This permitted us to analyze various peculiar situations. 

In the high temperature/large volume limit we immediately obtained the expected thermody- 
namic result: entropy is extensive and, for a massless field, goes as VT d , where V is the volume 
of the region and d the number of spatial dimensions. At leading order our results match the 
standard calculation of the entropy of a field system with appropriate conditions at its boundary 
[3]. In other words, in the NW approach, a generic subsystem of a larger thermal system is also 
approximately thermal and has the same temperature. Again, this seems to suggest that NW 
localization is more appropriate for the spatial coarse-graining of microscopic quantities. 
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At low temperature/small volume, thermal entropy in our calculation goes to zero, but, differ- 
ently from the standard approach, no unusual subtraction is needed. In the canonical approach 
[131 E] thermal entropy is sub-extensive, Stherm — (VT d )( d+1 ^ d at low temperature, whereas our 
regularized entropy approaches extensivity from above (S ~ — VT d hiVT d for small VT d ). Since 
the calculations in the two approaches follow two completely different routes, it is difficult to 
recognize the technical reason for this discrepancy. Note, however, that at very low temperatures 
the modes that are typically excited have wavelengths much larger than the size of the subsystem 
itself. Such modes are not contemplated in the spectrum of the locally defined Hamiltonian, whose 
lowest non-zero eigenvalues are of order ~ 1/V 1 ' d . Those long-wavelength correlations between 
the internal and the external region dominate in this limit and appear to be the cause of our 
super-extensive entropy. Note also that — elne is the generic asymptotic behavior of the Von 
Neumann entropy of a density matrix with a parameter e which becomes a pure state in the limit 
e 0. 

Our renormalization procedure, although unconventional, looks encouraging. The thermody- 
namical description of quantum field systems is a lively subject, which finds its application beyond 
the physics of complex systems, going through a better understanding of interconnections between 
gravitational and quantum effects, as for example black hole thermodynamics or AdS/CFT cor- 
respondence. However, many of such interesting applications are affected by the problem of the 
ultraviolet divergences so that, apart from some exceptional cases, many results remain at a qual- 
itative or conjectural level. If holographic entropy bounds (e.g. [27J[28]) and the area-dependent 
black hole entropy (see [29], and also [30, 31] for entropy renormalization in that context) are 
to be taken as meaningful signals of quantum gravity effects, one may want to consistently get 
rid of the comparable area dependent contribution that appears already at low energies in flat 
space. Because of the automatic regularization, the finiteness of our approach seems to provide a 
powerful and concrete computational method to overcome such technical empasses. 
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A Tracks of a proof 

Let us consider the integral 

t f a [a m m 1 ---m n+1 

J = L dx °- L " +1 K + (* - *oi • ■ ■ K +1 + - ' (84) 

where the ranges i?j can be the interval It, = [-L/ ft, L//3] or the set £"1, = (— 00, —L/ ft]U[L/ ft, 00), 
and n > 1 (the case n = can be verified by a direct computation). Without loss of generality 
we can assume R n +\ = El- We can perform the integration in dx n+ \ by means of the formula 

ab 



dx 



[a 2 + (x - y) 2 ][b 2 + (z-x) 2 } 

b[b 2 — a 2 + (z — y) 2 ] arctan 5=2 + a[a 2 — b 2 + (z — y) 2 ] arctan • 

= [(y - z) 2 + (a +l) 2 }[(y - z) 2 + (a - b) 2 } 



so that 



[(y - z) 2 + (a + b) 2 ][(y - z) 2 + (a - b) 2 ] 



, , / j m n m n+l 

W{m n ,m n+1 ;x n ,x ;ft) = / dx„ +r 



55) 



[m 2 + (x n - x n+1 ) 2 ][m 2 l+1 + (x n+1 - x ) 



L 

m n+1 [m 2 n+l -m 2 n + (x - x n ) 2 } arctan 
[(x n - x ) 2 + (m n + m n+1 ) 2 ][(x n - xq) 2 + (ra n - m n+ i) 2 ] 



"^n[^n — m n+i + ( x o ~ ^n) 2 ] arctan 



m n+ i 



[(x n - x ) 2 + (m n + m n+ i) 2 ][(x n - x ) 2 + (m n - m n+1 ) 2 } 
m n m n+ i(x n - x ) log : 



^+i+(#-^o) 2 

[{x n - x ) 2 + (m n + m n+l ) 2 }[(x n - x ) 2 + (m n - m n+1 ) 2 } 

-(Q _> _/3) _| 7r(m n+1 + m„,)[(m n+1 - m n ) 2 + (x - x n ) 2 } 

[{x n - x ) 2 + (m n + m n+1 ) 2 ] [(x„ - x ) 2 + (m n - m n+1 ) 2 ] ' 

Next, a very careful analysis is needed, distinguishing the cases if R , R n are of type II and/or 
El- Note that it is an odd function of ft, apart from a term which do not contain ft. Indeed such 
term imply W(m n , m n+ i; x n , Xq; ft) — > if ft — > + . A very lengthy computation shows that one 
can find a set of positive constants such that 

W(a, b; x n , x ; ft) < f (87) 

After substitution in (1811) . we see that in the worst casd^ </ takes the form (1501) so that the limit 
/3 — > + exists. 

7 that is when all the remaining insertions are the black ones 
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